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We calculate the double transverse-spin asymmetries, &&tt(Qt)i in transversely 
polarized Drell-Yan process at the transverse-momentum Qt of the produced lepton 
pair. We perform all-order resummation of the logarithmically enhanced contributions 
in the relevant Drell-Yan cross sections at small Qt, which are due to multiple soft 
gluon emission in QCD, up to the next-to-leading logarithmic accuracy. The asym- 
metries £#tt(Qt) to be observed in polarized experiments at RHIC and J-PARC are 
studied numerically as a function of Qt- We show that the effects of the soft gluon re- 
summation to the polarized and unpolarized cross sections largely cancel in £&tt(Qt), 
but the significant corrections still remain and are crucial for making a reliable QCD 
prediction of ^tt(Qt)- In particular, the soft gluon corrections enhance £^tt(Qt) 
considerably in the small Qt region compared with the asymmetry in the fixed-order 
a s perturbation theory. We also derive a novel asymptotic formula which embodies 
those remarkable features of £&tt(Qt) at small Qt in a compact analytic form and 
is useful to extract the transversity 5q(x) from the experimental data. 
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1 Introduction 



In recent years, a variety of experiments has been devoted to explore the spin-dependent 
phenomena in hard processes. Especially, experiments with transvesely polarized hadrons 
have opened a new window to study rich structure of perturbative/nonperturbative dynam- 
ics of QCD associated with the transverse spin [I]. One of the fundamental quantities which 
newly enter into play is the chiral-odd, twist- 2 parton distribution, called the transversity 
5q(x); it represents the distribution of transversly polarized quark inside transversly po- 
larized nucleon, i.e, the partonic structure of the nucleon which is complementary to that 
associated with the other twist-2 distributions, such as the familiar density and helicity dis- 
tributions q(x) and Aq(x). However, Sq(x) has not been well-known so far. This is because 
8q(x) cannot be measured in inclusive DIS in contrast to q(x) and Aq(x); the chiral-odd 
nature requires a chirality flip, so that Sq(x) must be always accompanied with another 
chiral-odd function in physical observables. It is very recent that the first global fit of 5q(x) 
is given [2] using the semi- inclusive DIS data, in combination with the e + e~ data for the 
associated chiral-odd (Collins) fragmentation function. 

Transversely polarized Drell-Yan (tDY) process, p*p* — > X, is another promising 
process to access the transversity 5q(x). Based on QCD factorization, the spin-dependent 
cross section A T da = (da**—da*^)/2 is given as a convolution, A T da = J dx\dx 2 8H(x\, x 2 , n 2 F ) 
A T da(x° l / X\, x\jx 2 , Q 2 , (i 2 F /Q 2 ), where Q is the dilepton mass, fi F is the factorization scale, 

SH(xx,x 2 ; /4) = ^2 e ? W x i> f^F^Qi^, I*f) + Sq(x u /i 2 F )5q(x 2 , fi%)] , (1) 

is the product of transversity distributions of the two nucleons, summed over the mass- 
less quark flavors q with their charge squared e 2 , and A F da = (da** — da*')/2 is the 
corresponding partonic cross section. x\ = yfr e y ',x 2 = \fr e~ y are the relevant scaling 
variables, where r = Q 2 /S, and \fS and y are the total energy and dilepton's rapidity 
in the nucleon-nucleon CM system. At the leading twist level, the gluon does not con- 
tribute to the transversely polarized, chiral-odd process, corresponding to helicity-flip by 
one unit. The unpolarized cross section, da = (da** + da*^)/2, obeys factorization sim- 
ilar as At da, in terms of H(x\, x 2 , fj, F ) that is given by flTJ with 5q — > q and Sq — > q, 
and additional functions involving the gluon distribution that comes in as higer-order a s 
corrections. Therefore, the double-spin asymmetry in tDY, Att = ATda /da, in principle 
provides clean information on the transversity 5q(x). At the leading order (LO) in QCD 
perturbation theory, x\ 2 coincide with the momentum fractions carried by the incident 
partons, e.g., ATda oc 5(x± — x\)5(x 2 — x 2 ), so that [3], [1] 

A T da 1 /n ^ SH(x° 1 ,x° 2 ;Q 2 ) + --- 
= da = 2 C ° S(20) H(x°,xW) + --- ' (2) 

where <j) denotes the azimuthal angle of one of the leptons with respect to the incoming 
nucleon's spin axis, and the ellipses stand for the QCD corrections of NLO or higher. The 
cos(20) dependence is characteristic of the spin-dependent cross section A^dcr of tDY [3]. 
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Att to be observed in tDY at RHIC-Spin experiment was calculated by Martin et al. [I] 
including the NLO QCD corrections. The results are somewhat discouraging in that the 
corresponding Att are at most a few percent [3].Q The reason is twofold (see (jSJ)): (i) tDY 
in pp collisions probes the product of the quark transversity-distribution and the antiquark 
one as (jTj), and the latter is likely to be small; (ii) the rapid growth of the unpolarized 
sea-quark distributions in H(x\,x\;Q 2 ) is caused by the DGLAP evolution in the low-x 
region that is typically probed at RHIC, VS = 200 GeV, Q < 10 GeV, and y/f < 0.05. 
Thus, small Att at RHIC appears to be rather general conclusion (see also [5]). 

We note that those previous NLO studies of A T t of fl2]) correspond to tDY with the 
transverse-momentum Qt of the produced lepton pair unobserved, and use the cross sec- 
tions ATda,da integrated over Qt in However, in view of the fact that most of the 
lepton pairs are actually produced at small Qt in experiment, it is important to examine 
the double transverse-spin asymmetries at a measured Qt, in particular its behavior for 
small Qt- This is defined similarly as (j2j) using the "QT-differential" cross sections, and we 
denote it as s^tt{Qt) distinguishing from the conventional Qr-independent Att- In fact, 
participation of the new scale Qt{^ Q) causes profound modifications of the relevant theo- 
retical framework. For example, now the numerator and the denominator of &/tt(Qt) ma y 
involve the parton distributions associated with the scales ~ Qt, such as 5H(x\, Q T ) 
and H(x1, x®; Q T ), respectively. For H{x\, x®; Q T ), the low-x rise of the unpolarized sea- 
quark distributions, mentioned in (ii) above, is milder compared with H(x\,x\; Q 2 ). Thus, 
if the former components play dominant roles compared with the latter in the denominator 
of &&tt{Qt) by certain partonic mechanism, &&tt{Qt) for small Qt region at RHIC can be 
larger than Att', the necessary partonic mechanism is indeed provided as the large loga- 
rithmic contributions of the type ln(Q 2 /Q T ), which is another remarkable consequence of 
the new scale Qt'- the small transverse-momentum Qt of the final lepton pair is provided 
by the recoil from the emission of soft gluons which produces the large terms behaving 
as a n s ln m (Q 2 / Q T ) / Q% (jn = 0, 1, . . . , 2n — 1) at each order of perturbation theory for 
the tDY cross sections. Actually, such enhanced "recoil logarithms" spoil the fixed-order 
perturbation theory, and have to be resummed to all orders in a s to make a reliable pre- 
diction of the cross sections at small Qt- Recently, we have worked out the corresponding 
"QT-resummation" for the tDY cross sections up to next-to-leading logarithmic (NLL) 
accuracy, which corresponds to summing up exactly the first three towers of logarithms, 
<ln m (Q 2 IQ 2 t)IQ 2 t with m = 2n - 1, 2n - 2 and 2n - 3, for all n Utilizing this result, 
in the present paper, we develop QCD prediction for £&tt{Qt) as a function of Qt- 

We will demonstrate that the soft gluon corrections are significant so that s^tt{Qt) in 
the small Qt region is considerably large compared with the known value for Att- HI I n 
addition to £^tt{Qt) in tDY at RHIC, we calculate £&tt{Qt) to be observed at J-PARC 
when the polarized beam is realized [7]. The latter case is also interesting because the 
fixed target experiments at J-PARC probe the parton distributions in the medium x region 

*In [4], the corresponding asymmetries are defined through certain integration over 0, and equal ([2|) 
with the formal replacement cos(2</>) — * 2/tt. 

^For the impact of the Qt resummation on the spin asymmetries in semi-inclusive deep inelastic scat- 
tering, see [8]. 
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(y/S = 10 GeV, Q>2 GeV, and a/t>0.2), and thus large asymmetries are expected even 
for the Qr-independent A TT j9] (see (ii) above). We also find that £#tt(Qt) f° r Qt ~ 
deserves special attention from theoretical as well as experimental point of view, and derive 
a compact analytic formula for s^tt{Qt ~ 0). 

The paper is organized as follows. In Sec. 2, the QT-resummation formula for the tDY 
cross sections is introduced, and all ingredients necessary for calculating the Q^-dependent 
asymmetries £#tt(Qt) including the NLL resummation contributions are explained. In 
Sec. 3, numerical results of s^tt{Qt) at RHIC and J-PARC are presented. Sec. 4 is devoted 
to the discussion of analytic formula of s^tt{Qt) at Qt ~ using the saddle-point method. 
Conclusions are given in Sec. 5. 



2 Resummed cross section and asymmetry for tDY 

Throughout the paper we employ the MS factorization and renormalization scheme with 
the corresponding scales, fip and /xr. We first recall basic points of the fixed-order calcu- 
lation of the spin-dependent, QT-differential cross sections of tDY [6]. In the lowest-order 
approximation via the Drell-Yan mechanism, the lepton pair is produced with vanishing 
Qt, so that the corresponding partonic cross section is proportional to 5(Q T ). The one-loop 
corrections to the partonic cross section involve the virtual gluon corrections, and the real 
gluon emission contributions, q + q — » l + l + g; in the latter case, the finite Qt of the lepton 
pair is provided by the recoil from the gluon radiation. Those have been calculated in the 
dimensional regularization [6], and the differential cross section of tDY is obtained as 

t&mm = cos W) 3ot [AtX (Qi ' Q2 ' y)+ AtY {Ql • Q2 • y)] • (3) 

where AtX and AtY are, respectively, expressed as the convolution of ([1]) with the 
corresponding partonic cross sections, see [6] for their explicit form in the MS scheme: 
A T X = A T X^ + A T X^ as the sum of O(a s ) and 0(a\) contributions, where a s = a s (//^) 
with fifj the renormalization scale, and AtX^ = 5H{x\, x° ; fJ/p) 5(Q T ). The partonic 
cross section associated with AtX^ contains all terms that are singular as Qt — > 0, be- 
having Q T 2 x (\ii(Q 2 /Q T ) or 1) or 5(Q T ), while the 0(a s ) terms that are less singular than 
those in AtX^ are included in the "finite" part AtY. In ([3]), AtX becomes very large as 
~ a s ln(Q 2 /Q T )/Q T and ~ a s /Q T when Qt Q, representing the recoil effects from the 
emission of the soft and/or collinear gluon, and those terms have to be combined with the 
large contributions of similar nature that appear in each order of perturbation theory as 
a^\n 2n ~ 1 (Q 2 /Q T )/Q T , <ln 2n ~ 2 (g 2 /Q 2 T )/Q 2 T , and so on, from the multiple gluon emission. 

The resummation of those logarithmically enhanced contributions to all orders has been 
worked out [6], in order to obtain a well-defined, finite prediction for the cross section. This 
is carried out by exponentiating the soft gluon effects in the impact parameter b space, up 
to the NLL accuracy. As the result, AtX of ([3]) is replaced by the corresponding NLL 
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resummed component as AtX — > AtX nll , with [6] 



POO 

A T X^ h {Q^Q\y) = V> 2 / db-J (bQ 

i,j,k Jo 1 



"0 



b 2 



,0 "0 



b 2 



(4) 

Here J^bQx) is a Bessel function for the two-dimensional Fourier transformation from 
the b space to the Qt space, and bo = 2e~ 7B with 7^ the Euler constant. The symbol 
(g) denotes convolution as (C^ ® f 3 ) (x,fi 2 ) = J (dz/z) Cij(z,a s ([i 2 )) fj(x/ z, (i 2 ). Note 
that the suffix i, j, k can be either q, q including the flavor degrees of freedom, and we set 
f q (x,fi 2 ) = Sq(x, /i 2 ), fq(x,fi 2 ) = 5q(x,fi 2 ). The soft gluon effects are resummed into the 
Sudakov factor e s( - b '^ with 



S(b,Q) 



« 2 dK 2 



2 K z 



A q (a s (K 2 ))\n% + B q (a s (K 2 )) 



H 



(5) 



The functions A q , B q as well as the coefficient functions Cij are perturbatively calculable: 



A q (a s 



\poo (a„\ n A( n ) D 



a, 



E~=i (§ )" 4 n) , and Cyfra.) = 6^6(1 - z) + 



E:=i (t) n 4 At the NLL accuracy, 
A« = 2CV, A^ = 2C f ' 



7T 

IT 



G 



G 



9' 



-N 



BP = -3C F 



(6) 



where Cp = (N 2 — l)/(2iV c ), C*g = iV c , and Nf is the number of QCD massless flavors, and 

r 2 



C ( i )\z)=5 l3 C F 



7T 



y-4H(l 



(7) 



are derived in [6]. The result (jfJJ) coincides with that obtained for other processes [TUl [TT] . 
demonstrating that { Aq 1 ^ , A q 2 \ Bq 1 ^ } are universal (process-independent). Substituting 
([6]) and the running coupling constant a s (K 2 ) at two-loop level, the k 2 integral in dSJ) can be 
performed explicitly to the NLL accuracy, and the result can be systematically organized 
as (see also |13l IT] 



S(b,Q) 



-h {0 \\) + /i (1) (A) 



where the first and second terms collect the LL and NLL contributions, respectively, as 



(8) 



A (i) 
2ir[3 2 

2tt/3 3 
1 



[A + ln(l-A)], 



(9) 



ln 2 (l - A) 



A + ln(l - A) 



4tt 2 /3 2 _ 



1 - A 

2 



Dq 



ln(l - A) 



4 2 )-2vr/Vl«ln% 



/4 



A 



A 



2vr/3 
+ ln(l - A) 



(10) 



*B q n) (n > 2) and C^^z) (n > 1) depend on the process [12]. Also, A\ n} (n = 1,2,.. .) and B ( q L> are 

independent of the factorization scheme, but B q n ^ (n > 2) and C^j(z) (n > 1) depend on the factorization 
scheme (see e.g. [14]). 
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In these equations, (3q , (3\ are the first two coefficients of the QCD (3 function given by 
A) = (HC G - 2JV / )/(12tt), fix = (17C 2 - 5C G N f - 3C F N f )/(2Air 2 ), and 

Q 2 b 2 

A — fi a s (iJL%) In —ro- = /3 a s (ii 2 R )L . (11) 

In the b space, L = \n(Q 2 b 2 /6g) plays the role of the large logarithmic expansion parameter 
with b ~ 1/Qt, and A of (fTTj) is formally considered as being of order unity in the resummed 
logarithmic expansion to the NLL in (151), where the neglected NNLL corrections are down 
by a s {n 2 R ). Note that, expanding the above NLL formula (jl]) with (loT)- (fTTl) in powers of 
a s (fx 2 R ), the first three towers of logarithms, a™ \n m (Q 2 / Q 2 ^) / Q?p with m = 2n — l,2n — 2 
and 2n — 3, in the tDY differential cross section are fully reproduced for all n. Combining 
this expansion with the finite part Aj-Y of ([3]), the result gives the tDY differential cross 
section which is exact up to 0(a s ); thus we use the NLO parton distributions in the MS 
scheme for fj(x,fi 2 ) in (J2J), as well as for those involved in A^Y. 

We explain some further manipulations for our NLL formula; those were actually per- 
formed in [6] , but were not described in detail. The integrand of (jl]) depends on the parton 
distributions at the scale bo/b, according to the general formulation |15j . Taking the Mellin 
moments of AtY nll (Q^, Q 2 , y) with respect to the DY scaling variables x\ 2 at fixed Q, 

^tX n n ^ 2 (Q 2 t ,Q 2 ) ee fdxl (x?)^" 1 \\x\ (4) N2 ^A T X^(Q 2 T ,Q 2 ,y), (12) 

Jo Jo 

the 6-dependence of those parton distributions can be disentangled because the moments, 

(x, /i 2 ), obey the renormalization group (RG) evolution as fi^ipl/b 2 ) 
= J2j Uij t N{bl/b 2 , Q 2 )fj,N{Q 2 ), where Uij^lfi 2 , fi' 2 ) are the NLO evolution operators for the 
transversity distributions which are expressed in terms of the corresponding LO and NLO 
anomalous dimensions [HI [T7] and the two- loop running coupling constant. For (|T2|) with 
(jl]) and the above RG evolution substituted, several "reorganization" of the relevant large- 
logarithmic expansion is necessary for its consistent evaluation over the entire range of Qt, 
following the systematic procedure in [T3] elaborated for unpolarized hadron collisions: ex- 
ploiting the RG invariance, we have C ij>N (a s (b 2 /b 2 )) = C ljiN {a s {Q 2 ))e [asi ^ )c ^ /2w]x/{1 - x) 
to the corrections down by a s (fi 2 R ), for the iV-th moment of the coefficient function of ([7]), so 
that we make the replacement Cjj i Ar(a s (feo/^ 2 )) — ► Cij t Ki{a s {Q 2 )) = <%[l+(«s(<3 2 )CV/47r)(7r 2 — 
8)], up to the corrections of NNLL level for (I12p . Similarly, performing the large-logarithmic 
expansion for explicit formula of the NLO evolution operator Uij^(bl/b 2 , Q 2 ), we find 

U ihN (bl/b 2 ,Q 2 ) = dye^W R N (X) ee ^ P f N ln(l - A), (13) 

up to the corrections down by a s (fi 2 R ) which correspond to the NNLL terms when substi- 
tuted into (1T21. p. Here A T P qq ^ N = -2C F [ip(N+l) +7^-3/4] is the JV-th Mellin moment 
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of the LO DGLAP splitting function for the transversity. As a result, (fI2j) is expressed as 



1 + Z^Crtf - 8) 



(Q 2 ,Q 2 ), (14) 



^,(03-, 2 ) = / <J (6g T )e^) + ^W + ^W , (15) 
Jo * 

where 5H NltN2 (Q 2 ) is the double Mellin-moments of SH^x^x®; Q 2 ) of ([TJ), defined similarly 
as (fl2j) . The complete dependence on is included in the exponential factor e s<yb, ^ +RN ^ ( A ) +Rjv 2 ( A ) 
through L = ln^ 2 ^ 2 /^ 2 ,), so that all-order resummation of the large logarithms L and the 
associated 6-integral in (115j) are now accomplished at the partonic level. 

We also mention some other "reorganization" , which is explained in [6] and is necessary 
in order to treat properly too short and long distance involved in the h integration of ([To]) : 
firstly, to treat too short distance Qb <C 1, we make the replacement 

L^L = ln(Q 2 b 2 /b 2 + 1) , (16) 

in the definition ( TTT1) of A, following [T4j ; note that the integrand of ( TT5l) depends on the 
large- logarithmic expansion parameter only through A (see ([8])- ([T0]) . (TLB]) ). This replace- 
ment allows us to reduce the unjustified large logarithmic contributions for Qb <C 1, due to 
L 3> 1, as L — > and e s<yb ^ +RN ^ ^ +Rn 2 W — > l ; while L and L are equivalent to organize 
the soft gluon resummation at small Qt as L = L + (9(l/(<56) 2 ) for Qb ^> 1. Secondly, the 
functions (jHj) and (TTUl) in the Sudakov exponent (jEJ) are singular when A = floa s (fj, R )L — > 1, 
and this singular behavior is related to the presence of the Landau pole in the perturba- 
tive running coupling a s (n 2 ) in QCD. To properly define the b integration of (TT51) for the 
corresponding long-distance region, it is necessary to specify a prescription to deal with 
this singularity [13] : decomposing the Bessel function in (fT5l) into the two Hankel functions 
as Jo(bQ T ) = (H^ibQr) + H^ 2 \bQT))/2, we deform the fe-integration contour for these 
two terms into upper and lower half plane in the complex b space, respectively, and obtain 
the two convergent integrals as \b\ — > oo. The new contour C is taken as: from b = to 
b = b c on the real axis, followed by the two branches, b = b c + e ±l9 t with t G {0, oo} and 
< 9 < 7r/4; a constant b c is chosen as < b c < bi, where b = b^ gives the solution 
for A = 1. Note, this choice of contours is completely equivalent to the original contour, 
order-by-order in a s (/x|,), when the corresponding formulae are expanded in powers of a s . 
Therefore, this contour deformation prescription provides us with a (formally) consistent 
definition of finite 6-integral of (TT5]) within a perturbative framework. 

We now denote (TT4]) . with the replacement (fT6!) and the new contour C in (fT5|) . as 
AtX^^Q 2 ,, Q 2 ), and also denote the double inverse Mellin transform of A T X^ L N2 (Q^, Q 2 ), 
from (Ni,N 2 ) space to {x\,x%) space, as A T X NLL (Q^,Q 2 ,y). Defining (see ©) 

A T Y{Q%, Q 2 , y) = A T X(Q 2 T , Q 2 , y) + A T Y{Q 2 T , Q 2 , y) - A T X NLL (Q 2 , Q 2 , y) ^ , (17) 

where AyA^ NLL (Q 2 ", Q 2 , 2/)|fo denotes the terms resulting from the expansion of the re- 
summed expression up to the fixed-order a s (/^), we obtain the final form of our differential 
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cross section for tDY with the soft gluon resummation as [6] 



A T da 



dQ 2 dQ' T dyd(f) 



cos(20) 



a 



3N C SQ 2 



A T X^ L (Q 2 T ,Q 2 , y ) 



A T Y(Q 2 T ,Q 2 ,y) 



(18) 



From the derivation explained above, the expansion of this cross section in powers of a s (fi 2 R ) 
fully reproduces the first three towers of logarithms, a n s In" 1 (Q 2 j Q T ) / Q\ with m = 2n — 
l,2n — 2 and 2n — 3, associated with the soft-gluon emission for small Qt (<S Q), and 
also coincides exactly with the fixed-order result ([3]) to 0(a s ). Therefore, this formula 
( JT8|) avoids any double counting over the entire range of Qt- Note that A T Y(Q T , Q 2 , y) 
of (1171) corresponds to the "modified finite component" in our resummation framework: 
because the first and the third terms in the RHS of (1171) cancel with each other for Qt <C 
Q, A T Y(Q T ,Q 2 ,y) is less singular as Q T — > than Q^ 2 x (\n(Q 2 /Q' T ) or 1) or 5(Q T ), 
see the discussion below (j3J). Combined with A T X^ oc 5(Q T ), this also implies that 
A T Y(Q T ,Q 2 ,y) is of order a s (jJL 2 R ). In fact, (JED coincides exactly with A T Y(Q T ,Q 2 ,y) if 
( TT6l) is not performed. 

Because of this "regular" behavior of AtY(Q t , Q 2 , y) as Qt — > 0, we may consider ( jTTT) 
as the definition for the region where Qt > 0; in this case, the first two terms correspond 
to © for Q T > 0, i.e., 



A T da LO 



dQ 2 dQ?pdyd<j) 



cos (2(f>) 



a 



3N C SQ 2 



A T X^ (Q 2 t ,Q 2 , 



A T Y(Q 2 T ,Q 2 ,y) , (19) 



which gives the formula for the LO QCD prediction of tDY at the large-Q^ region. There- 
fore, our formula (fT8j) is actually the NLL resummed part, with the contributions to 0(a s ) 
(the third term of ffT7|) ) subtracted, plus the LO cross section; we refer to ffTS]) as the 
"NLL+LO" prediction, which gives the well-defined tDY differential cross section in the 
MS scheme over the entire range of Qt- It is straightforward to see that the integral of (TL81) 
over Q T reproduces that of ([3]) exactly, because L = at b = (see also [H]). 

We can extend the above results to unpolarized DY by mostly trivial substitutions to 
switch from spin-dependent quantities to spin-averaged ones, e.g., by removing "A^" and 
making the replacement, 5H(x\, x 2 ; /i 2 ) — > H (x±, x 2 ; /i 2 ), cos(20)a 2 /(3iV c S'(5 2 ) — > 2a 2 /(3N c SQ' 2 
etc., in the above relevant formulae. The explicit form of the spin-averaged quantities, such 
as X(Q T , Q 2 , y), Y(Q T ,Q 2 ,y), as well as those corresponding to the coefficient functions 
Cij(z,a s ) in (j4j), can be obtained from the results in [T5l EEE]. A different point from the 
polarized case is that now the gluon distribution f g (x,fi 2 ) = g(x,[i 2 ) participates, so that 



the suffix i,j of the "spin-averaged C. 



can be 



u g" 



as well as 



<q,<f. 



This also im- 



plies that A T P qq ,N appearing in (fTBl has to be replaced by the Mellin moment of the LO 
DGLAP splitting functions for the unpolarized case, which involve the mixing of gluon, 
and the "new Uij^(b^/b 2 ,Q 2 ) n represent the corresponding "evolution matrix" that was 
discussed in [TH [T3]. On the other hand, the formulae (l8l- (ITDl) of the Sudakov exponent 
hold also for the unpolarized case, reflecting that the coefficients relevant at the NLL 
level are universal [TTI [T2l [131 [H]. We list explicit form of the relevant formulae for the 
unpolarized cross sections in Appendix. 
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Taking the ratio of (TTB"]) to the corresponding NLL+LO prediction for unpolarized differ- 
ential cross section, we obtain the double transverse-spin asymmetry in tDY, for transverse- 
momentum Qt, invariant-mass Q, and rapidity y of the produced lepton pair, and azimuthal 
angle <fi of one of the leptons, as 

*M<M = i cq 3W ) ^" LL «&, Q 2 >v) + <?,v) . (20) 

2 KVI X"^(Ql,Q\ V ) + YiQl,Q\y) 

To the fixed-order a s without the soft gluon resummation, ( 120]) reduces to the LO prediction 
of the asymmetry for Q T > 0, 

1 A T xm(Q* T ,Q*,y)\ 2 >0 + A T Y(Q*,Q 2 , v ) 

"4t(Qt) - - 2 c-W) xm{Ql , Q ,, y C i>0 + Y(Ql,Q>,y) • (21) 

as the ratio of (fHJj) to the corresponding unpolarized cross section. 



3 The asymmetries &/ T t{Qt) at RHIC and J-PARC 

We evaluate the asymmetries, derived in the last section, as a function of Qt- We use 
the similar parton distributions as in the previous NLO studies jl] of Qr-independent 
Att of (T5]): for the transversity Sq(x,Q 2 ) participating in the numerator of the asymme- 
tries, we use a model of the NLO transversity distributions, which obey the corresponding 
NLO DGLAP evolution equation and are assumed to saturate the Soffer bound [19] as 
5q(x,fil) = [q(x, [j,q) + Aq(x, Hq)]/2 at a low input scale /x — 0.6 GeV using the NLO 
GRV98 [2D] and GRSV2000 ("standard scenario") [2T] distributions q{x,/4) and Aq(x,fi%), 
respectively. The NLO GRV98 distributions q(x,Q 2 ), g(x,Q 2 ) are also used for calculating 
the unpolarized cross sections in the denominator of the asymmetries. 

It is known that the QT-spectrum of DY lepton pair is affected by another nonper- 
turbative effects, which become important for small Qt region [15]: we have obtained the 
well-defined tDY cross sections and asymmetries that are free from any singularities, with 
a consistent definition of the integration in ffTBI) over the whole b region. However, the in- 
tegrand of (lisp involving purely perturbative quantities is not accurate for extremely large 
|6| region in QCD, and the corresponding long-distance behavior has to be complemented 
by the relevant nonperturbative effects. Formally those nonperturbative effects play role to 
compensate the ambiguity that the prescription for the b integration in (TT5|) to avoid the 
singularity in the Sudakov exponent S(b,Q) of (18"|)- (ITD|) is actually not unique (see [T5]). 
Therefore, following [JSJ [THl E] , we make the replacement in ffTB]) as 

e S(b,Q) _^ e S(b,Q)-g NP b 2 ^ (22) 

with a nonperturbative parameter g^p- Because exactly the same Sudakov factor e s( - b '®' 
participates in the corresponding formula for the unpolarized noted above f[2"0j) . 

we perform the replacement (1221) with the same nonperturbative parameter g^p in the 
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Figure 1: The spin-dependent and spin-averaged differential cross sections for tDY: (a) 
ATda/dQ 2 dQTdyd<j) and (b) da/dQ 2 dQTdyd<j), as a function of Qt at RHIC kinematics, 
yfS = 200 GeV, Q = 5 GeV, y = 2 and = 0, with g NP = 0.5 GeV 2 . 

NLL+LO unpolarized differential cross section contributing to the denominator of (1217]) . 
This may be interpreted as assuming the same "intrinsic transverse momentum" of partons 
inside nucleon for both polarized and unpolarized cases, corresponding to the Gaussian 
smearing factor of (122]) . We use gNP — 0.5 GeV 2 , suggested by the study of the Qt- 
spectrum in unpolarized case [22]. 

For all the following numerical evaluations, we choose = for the azimuthal angle 
of one lepton, fiF = = Q for the factorization and renormalization scales and b c = 0, 
9 = g7j7r for the integration contour C explained below (TT6]) . 

First of all, we present the transvserse- momentum Qy-spectrum of the DY lepton pair for 
\/S = 200 GeV, Q = 5 GeV, and y = 2, which correspond to the detection of dileptons with 
the PHENIX detector at RHIC. The solid curve in Fig.[^a) shows the NLL+LO differential 
cross section ([18]) for tDY, multiplied by 2Qt, with g^p = 0.5 GeV 2 for (1221) . We also show 
the contribution from the NLL resummed component A<r^ NLL in (fTB"]) by the dot-dashed 
curve, and the LO result using (]T9]) by the dashed curve. Fig. [Jib) is same as Fig. [T](a) but 
for the unpolarized differential cross sections. The LO results become large and diverge as 
Qt — ► 0, while the NLL+LO results are finite and well-behaved over all regions of Qt- The 
soft gluon resummation gives dominant contribution around the peak of the solid curve, i.e., 
at intermediate Qt as well as small Qt- To demonstrate the resummation effects in detail, 
the two-dot-dashed curves in Figs. [U^a), (b) show the LL result which is obtained from the 
corresponding NLL result (dot-dashed curve) by omitting the contributions corresponding 
to the NLL level, i.e., /i (1) (A), R Nl (A) , R N2 (X) in ([15]) and a s (Q 2 )C F (n 2 - 8)/2tt in (ED for 
the polarized case (see (IE]) and the discussion below ([H])), and similarly for the unpolarized 
case. The LL contributions are sufficient for obtaining the finite cross section, causing 
considerable suppression in the small Qt region. On the other hand, it is remarkable that 
the contributions at the NLL level provide significant enhancement from the LL result, 
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Figure 2: The asymmetries £#tt(Qt) at RHIC kinematics, VS = 200 GeV, Q = 5 GeV, 
y = 2 and = 0: (a) &&tt{Qt) obtained from each curve in Fig. [TJ (b) The NLL+LO 
£&tt(Qt) of (120]) with (122]) using various values for Qnp- 

around the peak region for both polarized and unpolarized cases, and the effect is more 
pronounced for the former. Among the relevant NLL contributions, the "universal" term 
h^\\) produces similar (enhancement) effect for both (a) and (b) of Fig. (H while the other 
NLL contributions, associated with the evolution operators and the 0(a s (Q 2 )) coefficient 
functions (see e.g. (TT3l) . (1T4]) ). give different effects to the polarized and unpolarized cases. 

Fig. EJ^a) shows the double transverse-spin asymmetries in the small Qt region for tDY 
at RHIC, obtained as the ratio of the results in Fig. Ufa) to the corresponding results in 
Fig. Hlb) for respective lines, so that the solid curve gives the NLL+LO result (1201) . the 
dot-dashed curve shows the NLL result, 

^/NLL/yn \ 1 In ,^T^ NLL (Qti Q 2 -> V) /oo\ 

s&tt (Qt) = - cost 20) — ~ , 23) 

TT W ^ 2 1 V) X NLL (Q 2 T ,Q 2 ,y) ' 1 ' 



and the dashed curve shows the LO result ([21]) . The NLL+LO result is almost flat for 
Qt — > as well as around the peak region of the NLL+LO cross section in Fig. [1] This flat 
behavior is dominated by the NLL resummed components, and reflects the fact that the soft 
gluon emission effects resummed into the Sudakov factor e s ^ b,c ^ with ([S]) are universal to the 
NLL accuracy between the numarator and denominator of (12"01 . Slight increase of the solid 
line for Qt — > is due to the terms oc \n(Q 2 /Q T ) contained in the "regular components" 
AtY and Y in (120]) (see (1171) ). but such weak singularities which show up only at very small 
Qt will be irrelevant for most practical purposes. The LO result, obtained as the ratio of 
the two LO curves divergent as Qt — > in Figs. Q](a) and (b), gives the finite asymmetry 
for Qt > 0, but it does not have the flat behavior, i.e., decreases for increasing Qt, and is 
much smaller than the NLL+LO result. On the other hand, we note that the LL result, 
retaining only the resummmed components corresponding to the LL level, is given by (see 
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flU, (US}) 

<^{Qt) = \ c°b(^) ^;|^ - ^ , (24) 

which is independent of Qt, because the Qr-dependent factor (fl5l) with S^fr, Q) H-i^Ar^A) + 
Rn 2 (\) — > /i^ ^(A)/a s (Q 2 ) is common for both polarized and unpolarized cases. Namely the 
LL resummation effects cancel exactly between the numerator and the denominator in the 
asymmetry (1241) . As indicated in ( 1241) . the resulting value shown by the two-dot-dashed 
curve in Fig. [2]^a) coincides with the Qp-independent asymmetry (T5]) up to the NLO QCD 
corrections; note that Att — 4.0% including the NLO corrections similarly as [I] (see Table 
[1] below). However, we recognize that the soft-gluon resummation contributions at the NLL 
level enhances the asymmetry at the small Qr-region significantly, compared with the LL 
or fixed-order result. This is caused by the enhancement of the cross sections in Fig. [1] dis- 
cussed above, due to the universal h^(X) term and the other spin-dependent contributions. 
In particular, the evolution operators like ([TBI in the latter contributions allow the partic- 
ipation of the parton distributions at the scale b /b ~ Q T , and the components associated 
with those parton distributions indeed play dominant roles due to the mechanism embodied 
by the Sudakov factor e s ( b '®> of ( fl5l) . Combined with the different x-dependence between 
the transversity and density distributions as noted in (ii) above, the resulting enhancement 
arises differently between (a) and (b) in Fig. [I], and enhances the asymmetry as in Fig. [2]^a). 

In Fig. [2(b) we show the NLL+LO asymmetries £?tt(Qt) of (fJUD, with ([22} using 
various values of g^p- Here the solid curve is the same as the solid curve in Fig. [2](a), using 
9np = 0.5 GeV 2 . The result demonstrates that our NLL+LO asymmetry in the relevant 
small-Qr region is almost independent of the value of g^p in the range g^p = 0.3-0.8 GeV 2 . 
Although, at RHIC kinematics, the QT-spectrum from the spin-dependent cross section 
( TIB"]) with (1221) receives a sizable smearing effect in the relevant small-Qr region [6], the 
corresponding ^jvp-dependence is canceled by the similar dependence of the unpolarized 
cross section in the asymmetry ( 1201) . In our framework, such cancellation of the gNP~ 
dependence between the numerator and the denominator of ( [201) is observed for all relevant 
kinematics of our interest at RHIC, and also at J-PARC discussed below. However, we 
mention that too small value of g^p is useless in practice: the Gaussian smearing factor 
of (1221) for giyp = 0.1 GeV 2 is insufficient to suppress sensitivity to the extremely large \b\ 
region in (TT51) . so that the b integration receives the "inaccurate" long-distance perturbative 
contributions considerably at small Qt, which lead to unstable numerical behavior for 
Qt~1 GeV. For all the following calculations, we use g^p = 0.5 GeV 2 . 

Fig.[3]shows the NLL+LO asymmetries sx? T t(Qt) of §M at RHIC kinematics, \/S = 200 
GeV and various values of the dilepton invariant mass Q, using y = 2 and y = for (a) and 
(b), respectively; the dashed curve in (a) is the same as the solid curve in Figs. E](a), (b). 
For all cases in Fig. [31 we observe the typical flat behavior of £^tt{Qt) m the small Qt 
region, similarly as Fig. [2J On the other hand, £^tt{Qt) increases for increasing Q, and the 
value in the flat region reaches about 10% for Q = 20 GeV in Fig. [3](a). Such dependence 
on Q is associated with the small-x behavior of the relevant parton distributions: smaller 
Q corresponds to smaller x\ 2 — e ±y Q/\fS, so that the small-x rise of the unpolarized 
sea-distributions enhances the denominator of f[2"0"|) . We obtain larger &?tt{Qt) for y — 2 



12 



15r 



(a) y=2 



£10- 



5- 



"0 0.5 



Q=2GeV 

Q=5GeV 

Q=8GeV 

Q=15GeV 

Q=20GeV 



1 1.5 2 

Q T (GeV) 



2.5 



15 



V 



[/I 

< 



0, 



0.5 



1 1 1 

(b)y=0 


1 i 1 i 

- Q=2GeV 


--- Q=5GeV 


Q=8GeV 




^-uuev 




Q=20GeV 




, i , i 



1 1.5 2 
Q T (GeV) 



2.5 



Figure 3: The NLL+LO ^tt{Qt) of (J2DH with (J22J using = 0.5 GeV 2 at RHIC 
kinematics, yfS = 200 GeV, = with y = 2 and y = for (a) and (b), respectively. 





Q = 2GeV 


Q = 5GeV 


Q = 8GeV 


Q = 15GeV 


Q = 20GeV 


y/S = 200GeV 


y = 2 


3.3% 


4.0% 


4.9% 


6.5% 


7.4% 


y = 


3.5% 


3.7% 


4.4% 


5.9% 


6.9% 


= 500GeV 


y = 2 


1.8% 


2.0% 


2.4% 


3.4% 


4.0% 


y = 


2.2% 


2.1% 


2.4% 


3.2% 


3.8% 



Table 1: The Qp-hrdependent asymmetry Att of ([2]) including the NLO QCD corrections 
at RHIC kinematics. 



compared with the y = case, but the y-dependence of £/tt(Qt) is not so strong for all Q. 
For comparison, we also evaluate the Qr-hidependent asymmetry Att of ([2]) including the 
NLO QCD corrections and with the same nonperturbative inputs as those used in Fig. [3j 
The results are shown in Table (TJ and these exhibit similar behavior with respect to the 
Q and y dependence as that in Fig. [3j Note that we reproduce the NLO value of A TT in 
Table[T]when we integrate respectively the numerator and the denominator of the NLL+LO 
asymmetry (1201) over Q T for each curve of Fig. [3] (see discussion below (fl9l ). But the NLO 
Att are smaller by about 20% than the corresponding values of the NLL+LO £$tt{Qt) in 
the "flat" region at small Qt- This enhancement of &&tt{Qt) compared with Att arises 
from the soft gluon resummation at the NLL level, as discussed in Fig. [2]^a) above. 

Fig. |4] is same as Fig. 3, but for another RHIC kinematics with v^S* = 500 GeV. General 
behavior for the Qt, Q and y dependence is similar as that in Fig. [31 Comparing the curves 
with the same values of Q, y between Figs. |3]and|4l &/tt(Qt) are smaller for higher energy 
= 500 GeV than those for = 200 GeV. This reflects the smaller x° lj2 = e ±y Q/V^ 
for larger \/S, and the corresponding enhancement of the denominator in ( 1201) . Similarly 
as Fig. El the NLL+LO ^tt(Qt) in the fiat region of Fig. 1 are larger by 20-30% than 
the corresponding NLO Att shown in Table [H It is generally true, regardless of the 
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Figure 4: Same as Fig. [3j but for = 500 GeV. 



specific kinematics or the detailed behavior of nonperturbative inputs, that the NLL+LO 
£0tt(Qt) °f ( 120]) in the flat region is considerably larger than the corresponding NLO 
Att, because this phenomenon is mainly governed by the partonic mechanism associated 
with the soft gluon resummation at the NLL level, as demonstrated in Figs. [I], [2j On the 
other hand, apparently the absolute magnitude of both s$tt{Qt) and Att is influenced 
by the detailed behavior of the input parton distributions, in particular, by their small-x 
behavior at RHIC. For example, if we change the input parton distributions, explained 
above (J22D, from the NLO GRV98 and GRSV2000 distributions into the NLO GRV94 [23] 
and GRSV96 [21] distributions, the NLO values of Att become smaller by 30-40% than the 
corresponding values in Table [D We note that the latter distributions are the ones used in 
the calculation of [I] , and the small-x behavior of the transversity distributions, resulting 
from 5q(x,fil) = [q(x, //q) + Aq(x, (J>q)]/2 at the input scale /i , is rather different between 
those two choices of the distributions, reflecting that the helicity distributions at small x 
are still poorly determined from experiments. U 

Next we discuss tDY foreseen at J-PARC. Fig. 0(a) shows the Qt spectrum of the 
produced lepton pair for J-PARC kinematics, v^S* = 10 GeV, Q = 2 GeV and y = 0. The 
curves show the spin-dependent differential cross sections, and have the same meaning as the 
corresponding curves in Fig. (U(a). The double transverse-spin asymmetries are obtained 
as the ratio of the results in Fig. 0(a) to the corresponding results for the unpolarized 
differential cross sections, as shown in Fig. 0(b). We see that the results at J-PARC obey 
the similar pattern as those at RHIC shown in Figs. (H [2j the flat behavior is observed 
for the NLL+LO £&tt{Qt) at Qt — > as well as around the peak region of the NLL+LO 
cross section, and this is dominated by the NLL resummed components. Also the soft- 
gluon resummation contributions at the NLL level enhances the asymmetry at the small 
QT-region significantly, compared with the LL of (1241) or the fixed-order LO result. As a 



§ We thank H. Yokoya and W. Vogelsang for clarifying this point. 
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Figure 5: The tDY at J-PARC kinematics, VS = 10 GeV, Q = 2 GeV, y = and <p = 0. (a) 
The spin-dependent differential cross section Arda /dQ 2 dQTdyd(p using g^p = 0.5 GeV 2 . 
(b) The asymmetries s^tt{Qt) obtained by using each curve in (a). 





Q = 2GeV 


Q = 2.5GeV 


Q = 3.5GeV 


VS = lOGeV 


y = 


12.8% 


12.9% 


12.5% 


y = 0.5 


13.9% 


14.8% 


15.9% 



Table 2: Same as Tabled] but for J-PARC kinematics. 



result, we get £^tt{Qt) — 15% as the NLL+LO prediction around the flat region, which 
should be compared with the corresponding prediction Att — 12.8% for (T5]) including the 
NLO corrections (see Table [2]). The reason why we obtain much larger values of £&tt{Qt), 
and also of Att, than the RHIC case is the larger x\ 2 = 0.2 probed at J-PARC, where 
the transversities are larger and the unpolarized sea distributions are smaller. Another 
difference compared with the RHIC case is that the contribution of the "regular component" 
A T Y of (|T7|) in Fig. [5](a), and the associated increase of the solid curve as Qt — > in 
Fig. 0(b), due to the terms oc \n(Q 2 /Q 2 T ) in A T Y and Y of (ED}, are more pronounced, but 
the latter effect shows up only for Qt < 0.5 GeV. 

In Fig. [6] we show the NLL+LO asymmetries £&tt{Qt) of (120]) at J-PARC kinematics, 
VS = 10 GeV and various values of Q, with y = and y = 0.5 for (a) and (b), respectively; 
the solid curve in (a) is the same as the solid curve in Fig. [5^b). We observe the flat 
behavior of £&tt(Qt) in the small Qt region, where £^tt(Qt) — 15-20% and these values are 
significantly larger than the corresponding results for the Qr-independent, NLO asymmetry 
A T t of ([2]), shown in Table [2j We note that the dependence of £&tt(Qt), as well as A T t, 
on Q is weak in contrast to the RHIC case; recall that the rather strong Q-dependence in 
Figs. [31 H] was induced mainly by the growth of the unpolarized sea-distributions for the 
small Xi2i probed at RHIC. 
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Figure 6: The NLL+LO ^tt{Qt) of (EQD with ([221) using g NP = 0.5 GeV 2 at J-PARC 
kinematics, v^S* = 10 GeV, = with y = and y = 0.5 for (a) and (b), respectively. 



4 The saddle point formula 

In Sec. 3, we have observed the universal flat behavior of NLL+LO s^tt{Qt) of ( [201) at 
small Qt-i including the region around the peak of each DY cross section in the numerator 
and denominator of (T20|) for both RHIC and J-PARC cases. We have also demonstrated in 
Figs. [2] and [5] that those flat behavior is driven by the dominant effects from soft gluon re- 
summation embodied by the NLL resummed components A^A" NLL and X NLL in ([2"0"|) . As a 
result, the values of &?tt{Qt) obtained in the "flat region" of the corresponding experimen- 
tal data may be compared, to a good accuracy, with (|23|) . Still, extraction of transversity 
distributions through such analysis should be a complicated task compared with the usual 
fixed-order analysis: in the flat region of £&tt{Qt), the b integration in the resummed part 
(Hj) (see also (TH1) . ( [T5l) ) mixes up the parton distributions numerically with very large per- 
turbative effects due to the Sudakov factor shown in Figs. [Hand [5], as well as with another 
nonperturbative effects associated with g^p of f[22[) . Thus in each of the numerator and 
the denominator of f[20[) . the information on the parton distributions is associated with a 
portion of the large numerical quantity whose major part would cancel in the ratio of f[20[) . 
and this fact would obscure the straightforward extraction of the transversity using the 
above formulae like (TH1) . f[T5"[) . in particular with respect to its accuracy. 

We are able to derive a simple analytic formula which allows more direct extraction of 
the transversity distributions from the experimental data in the flat region of £&tt{Qt) and 
also clarifies the accuracy of the resulting distributions. For this purpose, we first note that 
the extrapolation of £^tt{Qt) in the flat region to Qt = corresponds to the case without 
the (experimentally uninteresting) weak enhancement at very small Qt due to the terms 
oc \og(Q 2 /Q T ) in the regular components A?Y and Y of ([201 . so that the resulting value 
is very close to the Qt — > limit of f[23[) in both RHIC and J-PARC cases (see Figs. [2]- 
[6]). Namely, £/tt L (Qt = 0) ma y be considered to give a practical estimate of the data of 
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■^tt(Qt) in the flat region with a good accuracy. Then, at Qt = 0, the region |6| ~ 1/Aqcd 
becomes important for the b integration of the relevant resummation formula (I15p . Note 
that we can treat "safely" such long-distance region, corresponding to the boundary of 
perturbative and nonperturbative physics, owing to that the nonperturbative smearing 
( 122"]) suppresses the too long-distance region \b\ ;§> 1/Aq CD , and that the dependence on 
a specific choice of gNP cancels in the asymmetries ^tt L (Qt = 0) as demonstrated in 
Fig. [2](b) (see also pjj below). 

For simplicity in the presentation, we fix as [Ir = Q in the following. In the relevant 
region \b\ ~ 1/Aqcd, we have |L| ~ ln(Q 2 /A QCD ) ~ l/a s (Q 2 ), i.e., |A| ~ 1 (see (USD, dU))- 
Because all logarithms, L and ln(Q 2 /AQ CD ), are counted equally large for Q Aqcd, the 
resulting contributions to ([To]) are organized in terms of a single small parameter, a s (Q 2 ), 
but with a different classification of the contributions in the order of a s (Q 2 ) from the 
usual perturbation theory that can be used in an other region, < |6| < 1/Q: as discussed 
below (fTT|) . when A = 0(1), the NLL contributions in the Sudakov exponent © produce the 
(9(1) effects in the resummation formula (I15p . while the NNLL contributions could yield the 
corrections of 0(a s (Q 2 )). Therefore, when the region \b\ ~ 1/Aqcd is relevant as Qt — > 
and we neglect the NNLL contributions in (1151) . the other contributions that correspond 
to the same order of a s (Q 2 ) in (I14p should be neglected for a consistent treatment, as 
[1 + a s (Q 2 )C F (7t 2 - 8)/2tt] -> 1, so thatH 

A T X^ L N2 (Q 2 T = 0, Q 2 ) = 5H Nl , N2 (Q 2 )I Nl , N2 (Q 2 T = 0, Q 2 ) . (25) 

We note that the contributions to the NLL resummation formula for the unpolarized case 
can be classified similarly in the Qt — > limit; in particular, the present classification 
implies that the gluon distributions decouple for Qt — > by negelecting the 0(a s (Q 2 )) 
contributions of the corresponding coefficient functions Cjj (see ( 1ST]) in Appendix). 

It is also worth noting that this classification coincides with the "degree approxima- 
tion" discussed in [15] : in general, if one wants to evaluate the cross section for Qt ~ in 
an approximation where any corrections are suppressed by a factor of [ln(Q 2 /AQ CD )]~ ( ^ Ar+1 \ 
one needs a "degree N" approximation; i.e., for the perturbatively calculable functions in 
the general form of resummation formula (jl]) with (JSJ), one needs A q to order a^ +2 , B q to 
order af +1 , Cy to order af, and the f3 function to order af +2 . This indicates that the 
NLL accuracy for a resummation formula corresponds to the degree approximation when 
the region Qt ~ is considered. In particular, this implies that the 0(a s ) contribution 
in the coefficient function Cy should be neglected for Qt ~ 0; on the other hand, that 
contribution is necessary to ensure the NLL accuracy for Qt > Aqcd in the classification 
based on resummed perturbation theory of towers of logarithms, a" lia 2n ~ 1 (Q 2 /Q' T )/Q 2 ^, 
a n s \n 2n - 2 {Q 2 /Q 2 T )/Q 2 T , and a n s \n 2n ~ 3 (Q 2 / Q 2 T ) / Q 2 T [IB] [HE]. 

Now we evaluate the b integral of (fT5l) at Qt = according to the above classification. 
We have to use the exponentiated form in the integrand of (fT5l) without Taylor expansion, 

™ In principle, we should use this classification also for the numerical calculations presented in Sec. 3 
when Qt ~ 0. But we did not make the corresponding replecement for the coefficient functions Cij at 
Qt ~ in the calculations of Figs. [TJ[6] If we performed that replacement, the NLL+LO ([20]) as well as 
NLL ([23]) asymmetries at Qt ~ in those figures would increase by about 5%. 
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because the region |A| ~ 1 is relevant (see f[8])- ([T3]) ) [26] . This type of integrals can be 
evaluated with the saddle point method: we extend the saddle point evaluation applied to 
the LL resummation formula with g^p = [25], US] into the case of our NLL resumma- 
tion formula (fT5|) with nonzero Qnp- The corresponding extension is possible based on the 
present formalism that accomplishes resummation at the partonic level. We note that pre- 
vious saddle-point calculations consider the case with g^p = to avoid model dependence 
for the prediction of the cross sections, but the resultant saddle-point formula is applicable 
only to the production of extremely high-mass DY pair and is practically useless (see e.g. 
[21)] [T51 12"7]). We find that, with nonzero gNP-, we can obtain a new saddle-point formula 
applicable to the RHIC and J-PARC cases; also, although the behavior of the cross sections 
are influenced by specific value of gNP, the asymmetries are not, as already noted above. 

When Q is large enough, a s (Q 2 ) < 1, so that the b integral in ffT5]) with ([221 and Qt = 
is dominated by a saddle point determined by mainly the LL term in the exponent (jHJ) of 
the Sudakov form factor [25l [15]. In this case, the contributions to the b integration from 
too short <C l/Q) and long distance (\b\ ^> 1/Aqcd) along the integration contour C 
explained below f[T6"]) are exponentially suppressed: this allows us to give up the replacement 
([IE]) ; also we may neglect the integration along the two branches, b = b c + e ±l9 t with 
t G {0, oo}, in C, when b c is sufficiently large but is less than the position of the singularity 
in the Sudakov exponent, b^. In fact, we can check numerically that the relevant integrand 
has a nice saddle point well below (above 0) for the kinematics of our interest. Then, 
changing the integration variable to A, given by ([HI) , we get (see Q, ([TU[) . ([TBI) ) 

I Ni ,n 2 (Q 2 t = 0, Q 2 ) = An2 f° (n2 , t C dAe-fW+^W+^fAH^CA) (26) 

4:Q 2 (3 a s {Q 2 ) 

where A c = (3 a s (Q 2 ) ln{Q 2 b 2 Jb 2 ) (< 1), and 



A h^(\) g_xpb 2 

(3 a s (Q 2 ) a s (Q 2 ) Q 2 



C(A) = --^7^ " ^ + ^e«^> . (27) 



An important point is that the ratio, [h^(X) + Rn^X) + Rn 2 (\)]/((\), actually behaves 
as a quantity of the order of a s (Q 2 ) in the relevant region < A < A c of the integration 
in ( [26]) . even for nonzero g^p — 0.5 GeV 2 . The precise position of the saddle point in the 
integral of ([26J) is determined by the condition, -C'(A) + /i (1) '(A) + R Nl '(X) + Rn 2 'W = 0, 
and we express its solution as A = Xsp + AXsp where Xsp is the solution of C(X) = 0, i.e., 



A q Xsp 9Npb 



1 H - ^— = i^^e/WQ 2 ) 28 

27r/5 l-A SP Q 2 1 ) 

is satisfied, and AXsp = [h^'(Xsp) + Rni(Xsp) + Rn 2 '(Xsp)]/("(Xsp) denotes the shift of 
the saddle point at the NLL accuracy. Evaluating ([26]) around A = Xsp + AAsp, we get 



N2 (0,Q 2 ) = [ —3 JjL e -C(WW) ) e R Nl ^s P ) + R N2 (x SP )^ (29) 



4Q 2 f3 a s (Q 2 ) V C"(X SP ) 
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to the NLL accuracy. Here the contributions from the third or higher order terms in the 
Taylor expansion of the exponent in (126}) about the saddle point A = \$p + AAsp, as well 
as the other terms generated by the shift AAsp, are found to give the effects behaving as 
0{a s {Q 2 )), i.e, are of the same order as the NNLL corrections, and thus are neglected, 
similarly as in (jl4|) . (1251) . according to the classification of the contributions at Qp = 0. 
Substituting (129}) into (1251) and performing the double inverse Mellin transformation to the 
x®) space, the result is expressed as the factor in the parentheses of (129"]) . multiplied 
by (HJ with the scale, fi F -> b /b S p where b S p = {b Q / Q)e Xspl{Woaa{Q2) \ because e RN i {Xsp \ 
e Rjv 2 (Asp) j n (1291) can be identified with the NLO evolution operators from the scale Q to 
b /b S p, to the present accuracy (see (IT3"|) ). The saddle-point evaluation of the corresponding 
resummation formula for the unpolarized case can be performed similarly, and the result is 
given by the above result for the polarized case, with the replacement 5H(x\, x^'-, b\/b 2 SP ) — > 
H(xi,X2', b^/b 2 SP ). The common factor for both the polarized and unpolarized results, 
given by the contribution in the parentheses of (129]) . involves "very large perturbative 
effects" due to the Sudakov factor, and shows the well-known asymptotic behavior [25J, 
~ (A^ CD /Q 2 ) aln ( 1+1 / a ) with a = A$p /(2ir/3o), for Q > Aq CD ; but this factor cancels out for 
the asymmetry. As a result, we obtain the Qp — > limit of (1231) as 

(Qt - 0) - - cos(20) H{xlx ,_ hl/hlp) , (30) 

which is exact, up to the NNLL corrections corresponding to the 0(a s (Q 2 )) effects. This 
remarkably compact formula is reminiscent of stfppiQp) of (124}) that retains only the LL level 
resummation, or the Qp independent asymmetry of (|SJ), but is different in the scale of the 
parton distributions from those leading-order results. Namely, our result (130]) demonstrates: 
in the Qp = limit, the all-order soft-gluon-resummation effects on the asymmetry mostly 
cancel between the numerator and the denominator of (13"0"]) . but certain contributions at 
the NLL level survive the cancellation and are entirely absorbed into the unconventional 
scale bo/bsp for the relevant distribution functions. 

The new scale b /bsp is determined by solving (128]) numerically, substituting = 2Cp 
from ([6]) and input values for Q and g^p, but it is useful to consider its general behavior: 
the LHS of ( 128]) equals 1 at Xsp = 0, decreases as a concave function for increasing Xsp, 
and vanishes at X SP = + /(2tt(3 )} = 0.6; while the RHS is in general much smaller 
than 1 at Xsp = 0, increases as a convex function for increasing Asp, and is larger than 1 at 
Asp — 1. Thus the solution of ( 128]) corresponds to the case with LHS = RHS ~ 1/2, more 
or less independently of the specific value of Q and gNP, so that we get b /b S p ~ b y/2gNP- 
This result depends only mildly on the nonperturbative parameter gNP, and suggests that 
one may always use b /bsp — 1 GeV, for the cases of our interest where Q is of several GeV 
and g^p ~ 0.5 GeV 2 as in Figs. HE The actual numerical solution of (125]) justifies this 
simple consideration at the level of 20% accuracy. This fact will be particularly helpful in the 
first attempt to compare (130]) with the experimental data so as to extract the transversity 
distributions. 

Our saddle-point formula (130]) embodies the characteristic features of the NLL soft gluon 
resummation effects on the asymmetries ^tt(Qt), emphasized in Sec. 3. In particular, our 
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derivation of (IBOl) demonstrates clearly the mechanism, which makes the parton distribu- 
tions at the low scale ~ Qt play dominant roles, and leads to the "enhancement" of the 
dot-dashed curve in Figs. [2] and [5j As noted in the beginning of this section, (IBOl) may 
be directly compared with the experimental value of the asymmetries £&tt(Qt), observed 
around the peak of the Qt spectrum of the corresponding DY cross sections. But there 
is one caution for such application. As seen from the above derivation, the parton distri- 
butions appearing in (|30|) are the NLO distributions up to the corrections at the NNLL 
level; e.g., the transversity distributions appearing in the numerator of (|30|) is obtained 
by evolving the customary NLO transversity Sq(x, Q 2 ) at the scale Q, to the scale b /bsp 
using f[TH]) that is the NLO evolution operators up to the NNLL corrections. Therefore, the 
formula (1BU1) can be used in the region where NNLL corrections are small; we know that 
the NNLL corrections at Qt ~ correspond to 0(a s (Q 2 )) effects, and should be negli- 
gible in general. However, such straightforward estimate might fail at the edge region of 
the phase space, e.g., at the small x region: because the relevant evolution operators ([TBI) 
actually coincide with the leading contributions in the large-logarithmic expansion of the 
usual LO DGLAP evolution!! (|30|) would not accurate when the NLO corrections in the 
usual DGLAP evolution are large compared with the contributions of (fT3|) . Such situation 
would typically occur in the region with small x\ 2 , corresponding to the case with large 
\fS. In Tabled we compare j^tt L (Qt = 0) using the numerical b- integration ("NB"), 
obtained as the Qt — ► limit of the dot-dashed curve in Figs [3(a) andE^a), with those 
using the saddle-point formula (IBUl) . For the latter we use b /bsp obtained as the solution 
of ([28]) with qnp = 0.5 GeV 2 , and consider the two cases for the parton distributions par- 
ticipating in ([30]) : "SP-I" uses the parton distributions which are obtained by evolving the 
customary NLO distributions at the scale Q, to b /bsp using the NLO evolution operators 
up to the NNLL corrections like ([TBI) ; "SP-II" uses the customary NLO distributions at 
the scale b /bsp- Here the "customary NLO distributions" are constructed as described 
above (12"2"|) . First of all, the results for SP-I demonstrate the remarkable accuracy of our 
simple analytic formula (1301) for both RHIC and J-PARC, reproducing the results of NB 
to the 10% accuracy. E3 On the other hand, the results for SP-II indicate that the NNLL 
corrections are moderate for large y/S at RHIC, while those are expected to be small for 
small VS at J-PARC. 

We propose that our simple formula ([BO]) is applicable to the analysis of low-energy 
experiment at J-PARC in order to extract the NLO transversity distributions directly from 
the data. On the other hand, (jB0|) will not be so accurate for analyzing the data at RHIC, 
but will be still useful for obtaining the first estimate of the transversities. We emphasize 
that such (moderate) uncertainty in applying our formula (IBUj) to the RHIC case is not 
caused by the saddle-point evaluation, nor by considering the Qt — > limit, but rather is 

II This fact also suggests that one may use the fixed value, bo/bsp — 1 GeV, in ([30]) for all Q (and gNp) 
rather than solving |28|) numerically for each different input value of Q, <?jvp, because the sensitivity of the 
LO evolution on the small change of the scale is modest. 

**If we use the fixed value, bo/bsp = 1 GeV, for all cases, instead of the solution of (|28|). the results in 
SP-I change by at most 5%, for both RHIC and J-PARC kinematics. The corresponding change in SP-II 
is by less than 5% for J-PARC, and by about 10% (15%) for Q = 2-8 GeV (Q = 15-20 GeV) at RHIC. 
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Vs = 


200 GeV, y = 2 




Vs = 


= 10 GeV, 


y = 


Q 


2GeV 


5GeV 


8GeV 


15GeV 


20GeV 


2GeV 


2.5GeV 


3.5GeV 


SP-I 


4.3% 


5.4% 


6.6% 


8.7% 


9.8% 


14.1% 


14.5% 


14.8% 


SP-II 


7.3% 


8.7% 


9.8% 


11.8% 


12.7% 


14.7% 


14.8% 


14.2% 


NB 


3.8% 


4.9% 


6.1% 


8.2% 


9.4% 


13.4% 


14.0% 


14.9% 



Table 3: The Q T -> limit of ^ T LL (Q T ) of ([23]) for RHIC and J-PARC kinematics. SP-I 
and SP-II are the results of the saddle-point formula flHOl for g NP = 0.5 GeV 2 , using the 
evolution operators from Q to bo/bsp, to the NLL accuracy and to the customary NLO 
accuracy, respectively. NB is obtained from the dot-dashed curve in Figs[3^a) and[6](a). 

inherent in the general Qt resummation framework which, at the NLL level, implies the 
use of the evolution operators (fT3l with the LO DGLAP kernel; more accurate treatment of 
the small-x region of the parton distributions relevant to the RHIC case would require the 
resummation formula to the NNLL accuracy, where the NLO DGLAP kernel participates 
in the evolution operators (fl3"|) from Q to bo/b (see e.g. [14"]). 

5 Conclusions 

In this paper we have presented a study of double transverse-spin asymmetries for dilepton 
production at small Qt in pp collisions. The logarithmically enhanced contributions, which 
arise in the small Qt region due to multiple soft gluon emission in QCD, are resummed to 
all orders in a s up to the NLL accuracy. Based on this framework, we calculate numerically 
the spin-dependent and spin-averaged cross sections in tDY, and the corresponding asym- 
metries £^tt{Qt), as a function of Qt at RHIC kinematics as well as at J-PARC kinematics. 
The soft gluon resummation contributions make the cross sections finite and well-behaved 
over all regions of Qt-, so that the singular Qt spectra in the fixed-order perturbation theory 
are redistributed, forming a well-developed peak in the small Qt region. As a result, both 
the polarized and unpolarized cross sections become more "observable" around the pro- 
nounced "peak region" at small Qt-, involving the bulk of events. Reflecting the universal 
nature of the soft gluon effects, those large resummation-contributions mostly cancel in the 
cross section asymmetries £/tt{Qt), leading to the almost constant behavior of s^tt[Qt) 
in the small Qt region, but, remarkably, the effects surviving the cancellation raise the cor- 
responding constant value of £^tt{Qt) considerably compared with the asymmetries in the 
fixed-order perturbation theory. We have obtained a QCD prediction as £^tt{Qt) — 5-10% 
and 15-20% in the "flat region" for typical kinematics at RHIC and J-PARC, respectively, 
where the different values of s^tt{Qt) are associated with the different values of parton's 
momentum fraction probed by these two experiments. 

We have also derived a new saddle-point formula for s^tt{Qt ~ 0), clarifying the 
classification of the contributions involved in the resummation formula for Qt — > 0. The 
formula is exact to the NLL accuracy, and embodies the above remarkable features of 
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soft gluon resummation effects at small Qt in a compact analytic form. Our saddle-point 
formula may be compared with the data of £#tt{Qt) i n the peak region of the DY Qt- 
spectrum, and thus provides us with a new direct approach to extract the transversity 
distributions from experimental data. 

We mention that there is another kind of logarithmically enhanced soft-gluon contribu- 
tions, subject to the so-called "threshold resummation", besides those treated by the Qt 
resummation. It is known that the threshold resummation effects on the cross sections can 
be important when the probed momentum fractions of partons are rather large like at J- 
PARC. The corresponding effects for tDY are studied [28J in pp collisions at GSI kinematics, 
and the results indicate that the threshold resummation effects will not be so significant for 
the kinematical regions corresponding to experiments at J-PARC, and, furthermore, will 
cancel mostly in the asymmetries. 

We have revealed that the "amplification" of the double transverse-spin asymmetries 
£#tt(Qt) at small Qt is driven by the partonic mechanism participating at the NLL level, 
as the interplay between the large logarithmic gluon effects resummded into the universal 
Sudakov factor and the DGLAP evolutions specific for each channel. Thus similar phe- 
nomenon is anticipated also in pp collisions at the future experiments at GSI [29], where 
the large values are predicted for the Qt independent asymmetry Att [28j [30]. The appli- 
cation of our Qt resummation formalism to pp collision will be presented elsewhere |31j . 
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Appendix: Resummed cross section for unpolarized DY 

In this appendix, we summarize the corresponding results for the unpolarized Drell-Yan 
process which are necessary to calculate the asymmetry j&tt{Qt)- Although all of them 
have already appeared in the literatures jT8[ [15], we will list them for the convenience of 
the reader. 

The NLL resummed component, corresponding to (j3J) for the polarized case, reads, 
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where S(b, Q) is the exponent of the Sudakov factor, given by fl5])- ffTTj) . and 



7T 



C$(z) = C$>(z) = 2T R z(l-z) 



with Tr — 1/2. The parton distributions are defined as, 

H(xi , x 2 ] Hp) = ^2 e q [q( x i> Vf) Qfai Vf) + fai, t4) Q( x *> Pf)] 

q 

K x {x x , x 2 ; Hp) = ^2 e g M x i> Pf) + A*!-)] 9(%2, /4) » 



The finite component, corresponding to ApY of ([3]) for the polarized case, is given by, 

a s (/4) 



2tt 



[C F Y q (Q 2 T ,Q 2 ,y)+ T R Y g (Q 2 T , Q 2 , y 



where the first term comes from the qq annihilation subprocess as 



Y q (Q 2 T ,Q\y) 



2 

S 



dx\ H{x\ , x\\ Hp) 



v ^>e« x l x l 



X\X\ 



dx 2 H(x\ , x 2 ; Hf) 



-V X 2 Xn 



x\x 2 



+ 



Q 2 



T 



dxi Hi + I dx 2 H 2 + 2 H(x\ , x° 2 ] hD m 



(l-x+)(l-x+) 
(l-x°)(l-*°) 



with the variables according to [18j [6] O 



Q 2 + Q 2 T \ 1/2 
S J 



, } ■-) 



Q 2 + Q \ 
S 



2 \ 1/2 



X, 



x 2 x± — x\x 2 



Qt 



x 2 x 2 X\ x-i v S V S 

We used the shorthand notation for the integral that vanishes for Qp = 0, 



dxi Hi 



dxi 



~e« X l X l 
1 



H(x 1 , x* 2 ; Hf) < 1 + 



-2iy(x°,x°; a4) 



dxi 



.0 Xi — X{ 



H(xi , x 2 ; Hp) { 1 



x 



■I' i 



0\ 2 



2tf(x°,x°; //|) 



^Note that 2 = y / 7 r ^e ±y = a;* 2 = 2C? , 2 when = 0. 
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and J dx 2 H 2 = J dx\ H\ 
have 



y^-y 



Y g {Q 2 T ,Q\y) 



1 

S 



For the contribution from the gluon Compton process, we 



dx x T ^ , _„ 2 ^x x x 2 - t 



Kx(xi , x* 2 ; n F ) 

y/r+ey x l x \ 

1 



(xix^ 2 



Q 2 ! 



dx-) 



r+e 



-y X 2 X<i 



^/T+e-v x 2 x 2 

Ki(x\ , x 2 \ \i 2 F 



dX2 K^xx&f*** 



x\x 2 y 



( x 2 xf — r 2t(x 2 x~1 — t)' 



V 



x\x 2 



(x\x 2 y 



1 dx 2 

x° x 2 



+ 

+ (1 <-> 2) 



dxi 



t+ e 



y X\ X | 



K\{x\ , x 2 ; fj.%) |l — ~^~2 — ~| 
x 2 xf-r 2r(xlxf-r) 2 



X\X 2 



where 1 «-»• 2 means the exchange of the suffix of variables as well as y 
parton distributions, this exchange should be read, 



{xix* 2 f 



y. As for the 



K 1 (x 1 , xl) <-> K 2 (x\ , x 2 ) , Kx{x\ , x 2 ) <-> K 2 (xx , x* 2 ) , Jfi(x? , x 2 ) <-> ^2(^1 , ^ 
The fixed-order cross section in the MS scheme is given as 
da FO 2a 2 



dQ 2 dQ 2 T dyd<p 3N C SQ< 



[X (Q 2 T ,Q 2 ,y) + Y(Q 2 T ,Q 2 ,y)} 



(32) 



to the 0(a s ) accuracy, where X (Q 2 -, , Q 2 , y) = X NLL (Q 2 -, , Q 2 , j/)| FO denotes the "singu- 
lar" part resulting from the expansion of the resummed component up to the fixed-order 
a s (/x^). The formula (152]) should be compared with Q. Taking the similar steps as those 
in (fT2|) - (|T7|) . we obtain the differential cross section for unpolarized DY with the soft gluon 
resummation as 



da 



2a 2 



dQ 2 dQ 2 T dydct) 3N C SQ< 
which is used to calculate fl20l). 



X N ^(Q 2 T ,Q 2 ,y)+Y(Q 2 T ,Q 2 ,y) 



(33) 
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